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ABSTRACT 

Wc present self-consistent models for the radiative transfer in Shakiira-Sunyaev accre- 
tion discs in bright low-mass X-ray binaries (LMXB). Our calculations include the full 
effects of incoherent Compton scattering and the vertical temperature structure within 
the disc, as well as the effects of Doppler blurring and gravitational redshift. We find 
that the observed X-ray spectra are well fit by exponentially cutoff power-law models. 
The difference between the observed total spectrum and our calculated disc spectrum 
should reveal the spectrum of the disc/neutron star boundary layer and other emitting 
regions considered to be present in LMXB. 
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1 INTRODUCTION 



I The energy spectra of low-mass X-ray binaries (LMXB) are 
yf^ observed to be approximately exponential in shape with 
characteristic temperatures of a few keV. As the magnetic 
field of the central neutron star is weak, there is much emis- 
sion from the accretion disc, as well as from a boundary 
I layer between the disc and neutron star. Traditionally the 
O disc emission has been modelled by summing blackbody or 
modified blackbody spectra assuming that there is a uni- 
form temperature at each disc radius (Mitsuda et al. 1984; 
White et al. 1986; Czerny, Czerny & Grindlay 1986; Taam 
& Meszaros 1987; White, Stella & Parmar 1988; Nannurelli 
k, Stella 1989). However, such spectra must be poor ap- 
proximations to a real disc in which there is a vertical tem- 
perature gradient at each radius (and photons of different 
energies emerge from different depths) and in which the ef- 
fects of incoherent Compton scattering have been shown to 
be very important. 

Recenty, Shimura &; Takahara (1995) have presented 
self-consistent calculations for the radiative transfer and ver- 
tical structure in accretion discs of black hole candidates. 
They have extended the method that wc used (Ross, Fabian 
k, Mineshige 1992) for treating the radiative transfer in AGN 
accretion discs to include the vertical density structure as 
well as temperature structure within the disc. One of the 
models that Shimura & Takahara have presented was for a 
1.4-M0 central object, appropriate for a LMXB accretion 
disc. However, their total disc spectra were simple sums of 
locally emergent spectra; they have not been corrected for 
gravitational redshift and Doppler blurring to correspond to 
the spectra that would actually be observed. 



Here we present self-consistent models for the radiative 
transfer in Shakura-Sunyaev (1973) accretion discs about 
a 1.4-M0 neutron star, including first-order corrections for 
gravitational redshift and Doppler blurring. We make the 
standard assumption of a uniform density structure at each 
radius within the disc. Shimura & Takahara (1993) have 
found that, at high accretion rates, the density in the inner 
accretion disc is nearly uniform with depth (except in the 
outermost Thomson mean free path). We give simple fits to 
our final models in terms of exponentially cutoff power-law 
spectra. Wc hope that the fits will be useful for the interpre- 
tation of existing (e.g., Ginga) and future (e.g., XTE) data 
from LMXB. 

This work is part of an ongoing project to calculate self- 
consistent spectra of accretion discs. We plan to extend this 
work in the future to include the effects of heavy-element 
atomic processes and external illumination of the surface 
for LMXB accretion discs. 



2 METHOD 

We treat the radiative transfer in the inner, radiation- 
dominated portions of a LMXB accretion disc. A geomet- 
rically thin, Q-viscosit}' disc is assumed as described by 
Shakura &; Sunyaev (1973). The equations for the struc- 
ture of the disc have been summarized by Ross, Fabian 
& Mineshige (1992). The central neutron star has mass 
M = 1.4Mq. The accreting plasma is taken to consist of 
fully-ionized hydrogen and helium, with the relative abun- 
dance of helium given by Morrison &; McCammon (1983). 
Assuming that the accretion disc extends all the way down 
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to the last stable orbit at radius R = SRs = 6GM/c^ with- 
out striking the neutron star, the efficiency is 77 = 0.083, and 
the Eddington limit corresponds to an accretion rate 



Medd = 3.1 X 10" 



(1) 



These values are consistent with the Newtonian dynamics 
assumed by Shakura & Sunyaev (1973). 

At a given radius R, the density in the accretion disc 
is taken to be uniform in the vertical direction (see Laor &; 
Netzer 1989). Choosing dimensionless parameters 
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the flux emerging from the surface of the disc is given by 
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The density of the gas is given by 

3.4 X IQ-'r'^/- _3 
Po = ; -2 g cm , 
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and the half-thickness of the disc (measured in units of the 
Schwarzschild radius) is given by 



h 



(5) 



A value Q = 0.1 is assumed for the viscosity parameter. 

At each radius in the disc, the radiative transfer is 
treated using the Fokker-Planck/diffusion equation of Ross, 
Weaver &; McCray (1978) in plane-parallel geometry. 
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Here n is the photon occupation number, e is the photon 
energy, z is the vertical height within the disc, A^e is the 
free-electron number density, T is the gas temperature, js 
is the free-free spectral emissivity, k is the total opacity (per 
volume) and kb is the free-free absorption opacity. The first 
term on the right-hand side of equation (6) treats the full ef- 
fects of incoherent Compton scattering, including the lowest 
order Klein-Nishina correction to the scattering rate. The re- 
maining terms account for vertical diffusion of radiation as 
well as emission and absorption by free-free processes. The 
equation is solved in finite-difference form over a rectangular 
grid of e and z values as described by Ross, Weaver &; Mc- 
Cray (1978). The outer boundary condition for the radiation 
emerging from the surface of the disc has been discussed by 
Foster, Ross & Fabian (1986). 

The inner boundary condition employed in solving 
equation (6) deserves consideration. At a given radius, 
the radiative transfer could be treated for the entire half- 
thickness of the disc; the inner boundary condition would 
then consist of setting the spectral fiux equal to zero (F^ = 
0). In order to reduce the computing time, however, we only 
treat a portion of the disc thickness as described by Ross, 
Fabian k, Mineshige (1992). Care is taken that the "surface 
layer" treated is effectively thick, so that 



■\/3rff(TT -I- Tff) > 1, 



(7) 



even for high-energy photons. The spectral flux entering the 
surface layer from below is then assumed to satisfy the LTE 
radiative difi^usion equation (see Rybicki k, Lightman 1979), 

477 dB^{T) dT 



(8) 
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where Bi, (T) is the Planck function. The LTE temperature 
at the base of the surface layer is estimated as described 
by Ross, Fabian & Mineshige (1992), while the local tem- 
perature gradient is chosen so that the total flux entering 
the surface layer is consistent with the total rate of dynamic 
heating in the gas beneath the surface layer. In this way, the 
a priori uncertainty in the temperature at the base of the 
surface layer has very little effect on the calculated emergent 
spectrum. 

The occupation numbers, n(e, z), are relaxed until equa- 
tion (6) is satisfled throughout the grid of e and z values. 
During the relaxation, the local temperature, T{z), is found 
by balancing the total heating rate due to dynamic heating, 
Compton scattering, and free-free absorption with the total 
cooling rate due to inverse Compton scattering and free- 
free emission. Thus completely self-consistent, non-LTE so- 
lutions are found for the radiation fleld and the temperature 
structure of the gas. 



3 RESULTS 

It should be noted that the parameter / gives the fraction of 
the Eddington luminosity produced solely by the accretion 
disc. At least an equal luminosity must be produced in a thin 
boundary layer near the surface of the neutron star as the 
infalling gas loses its orbital kinetic energy (assuming that 
the neutron star is rotating slowly). Therefore, we limit our 
consideration to values / < 1/2. Sunyaev & Shakura (1986) 
have shown that the luminosity of the boundary layer can 
be as high as several times the luminosity of the accretion 
disc, depending on the actual radius of the neutron star. For 
/ < 1/2, the inner accretion disc also remains geometrically 
thin (h <^ r) as assumed by Shakura &; Sunyaev (1973). 



3.1 Local spectrum at r = 7 

First, we set / = 0.3 and consider the local spectrum emerg- 
ing from the accretion disc at r = 7. This is chosen as a char- 
acteristic spectrum since the quantity FoR^, which measures 
the contribution to the total disc luminosity, peaks at r « 7. 
For this model, the midplane of the disc is at a Thomson 
depth TT « 550. To assure that high-energy photons are 
thermalized at the base of the surface layer, we treat the 
radiative transfer down to a Thomson depth tt ~ 275. The 
calculated photon flux emerging from the surface of the disc 
is shown in Fig. 1. Also shown is a blackbody spectrum at 
the effective temperature Ten = (Fo/ct)^/"* = 8.2 x 10** K. At 
low photon energies (e ^ 0.2 kcV), the calculated flux mim- 
ics a blackbody at the surface temperature, which is some- 
what lower than Tes- At higher energies, however, the effects 
of Compton scattering become important. In the soft X-ray 
regime, the photon spectrum is nearly flat, while most of the 
flux emerges as a "Wien tail" in hard X-rays (e ^ 2keV). 
The hard X-ray tail results from photons escaping from 
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Figure 1. Calculated emergent spectrum (solid curve) at r = 7 
for / = 0.3. Shown for comparison are the elfective blackbody 
spectrum at Teg = 8.2 x 10^ K (dashed curve) and a Wien-law 
spectrum at Twien ~ 1-37 x lO'^ K with the same total energy flux 
(dotted curve). 

deeper, hotter layers of the accretion disc. The competition 
between Compton scattering of higher-energy photons and 
inverse Compton scattering of lower-energy photons results 
in a Wien-law spectrum, 

oc <syi^{-hvlkT^\^^). (9) 

To estimate the Wien temperature, Twien, wc examine the 
calculated temperature structure of the disc, seeking the 
Thomson depth for which the Compton y parameter equals 
unity (see Rybicki &; Lightman 1979): 

y = —tt = 1- (10) 

meC' 

The resulting temperature should represent the lowest for 
which Comptonization is complete and, therefore, is ap- 
propriate for the Wien law. For the model with / = 0.3 
and r = 7, equation (10) is satisfied at « 11 where 
T = 1.37 X lO'^K. Assuming that the entire energy flux 
emerges as a Wien spectrum with this temperature, we get 
the result shown in Fig. 1. It is seen to be in excellent agree- 
ment with the high-energy tail of the actual emergent spec- 
trum. Note that although high-energy photons are emitted 
at great depths where the temperature is higher, the Wien 
tail has a lower temperature corresponding to the smaller 
depth where y = 1. Without the effects of incoherent Comp- 
ton scattering, the emergent spectrum would be even harder. 
As it is, the Wien tail greatly exceeds the effective blackbody 
spectrum for e ^ 5 keV. 

The emergence of a Wien-law spectrum at high ac- 
cretion rates was predicted by Shakura & Sunyaev (1973). 
Comptonization should be complete or "saturated" when- 
ever the thermalization depth (where t» = 1) for high-energy 
photons is greater than the Comptonization depth (where 
y = 1). In terms of our parameters r and /, Shakura & 
Sunyaev predicted a Wien spectrum for 

r < 93a'/«/'^'- (11) 
For a = 0.1 and / = 0.3, this gives r ^ 25, which accounts 



Figure 2. Calculated emergent spectra at r = 7 for / = 0.1 
(lower solid curve), 0.2 (lower dashed curve), 0.3 (upper dashed 
curve), and 0.4 (upper solid curve). Also shown are a Wien-law 
spectrum (dash-dot curve) at Twien = 1-06 X 10'' K with the 
same total energy flux as the / = 0.1 model and the diluted 
blackbody spectrum (dotted curve) with q — 1.8 corresponding 
to the / = 0.4 model. 



for nearly 3/4 of the total luminosity of the disc. Only farther 
out in the accretion disc did Shakura & Sunyaev predict 
that Comptonization would be incomplete, resulting in a 
"modifled blackbody" spectrum. 

Calculated emergent spectra at r = 7 for models with 
accretion rates corresponding to / = 0.1, 0.2, 0.3, and 0.4 
are shown in Fig. 2. Even for / = 0.1, equation (11) predicts 
that Comptonization is complete for r ^ 12, accounting for 
half the total luminosity of the disc. For / = 0.1 and r = 
7, our calculations show that equation (10) is satisfied at 
Tr « 12 where T = 1.06 x lO'' K. A Wien spectrum at this 
temperature containing the entire flux is shown in Fig. 2, 
and it is seen to be in reasonable agreement with the high- 
energy tail of the emergent spectrum. 

A prominent Wien shoulder can be seen at e fa 3 keV for 
the highest accretion rate (/ = 0.4) shown in Fig. 2. Shimura 
&i Takahara (1995) used a "diluted blackbody" spectrum, 

= ^nB,{qTes), (12) 

to describe the Wien-law dominated portions of the local 
emergent spectra that they calculated. Here q is an appro- 
priate spectral hardening factor which, for high accretion 
rates, lies in the range 1.7 ^ q ^ 2.0. In Fig. 2 we show the 
diluted blackbody spectrum for the / = 0.4 model assuming 
a value q = 1.8. It is seen to be in good agreement with the 
calculated spectrum for e ^ 2 keV. 

Since we are interested in modelling high-luminosity 
LMXB accretion discs, we only consider high accretion rates 
(/ > 0.1). Accurate treatment of the Wien-law tail is im- 
portant in predicting the hard X-ray spectrum of such discs. 
Furthermore, the transition of the spectrum from a black- 
body in the EUV to a Wien tail in hard X-rays can only be 
modelled accurately by a fully self-consistent treatment of 
the radiative transfer, as we do here. 
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Figure 3. Summed spectra for an accretion disc witli / = 0.3. 
Upper solid curve sliows tlie total photon flux assuming a face- 
on disc at 1 kpc. Other curves show contributions to the total 
spectrum by different annuli assuming incoherent scattering for 
3 < r < 70 (solid curves), coherent scattering for 70 < r < 200 
(dotted curves), and blackbody spectra for 200 < r < 600 (dash- 
dot curve). Also shown is the total spectrum if the entire disc is 
assumed to emit blackbody spectra (dashed curve). 



3.2 Summed disc spectra 

We divide the region 3 < r < 200 of the accretion disc into 
20 annuli that make comparable contributions to the total 
luminosity. The typical emergent spectrum for each annu- 
lus is calculated and multiplied by the area of the annulus 
to find the contribution to the spectral luminosity. The dif- 
ferent contributions are added together, and the spectral 
photon flux is found assuming a face-on disc at an arbi- 
trary distance of 1 kpc. The result for an accretion disc with 
/ = 0.3 is shown in Fig. 3. For 17 annuli covering 3 < r < 70, 
a complete treatment of incoherent Compton scattering is 
performed as described previously. For three annuli cover- 
ing 70 < r < 200, however, y ^ 1 at the thermalization 
depth. Therefore, we treat Compton scattering as being co- 
herent for these calculations by dropping the Fokker-Planck 
term on the right-hand side of equation (6). Each coherent- 
scattering model is iterated until a self-consistent tempera- 
ture is found for use in the inner boundary condition (equa- 
tion 8). Finally, to give the spectrum a realistic behavior 
at very low photon energies, the spectrum for the region 
200 < r < 600 is found assuming that the disc emits an 
effective blackbody spectrum locally. 

The resulting calculated disc spectrum should be very 
accurate for e ^ 0.6 keV. The dip at e a 0.4 keV is probably 
somewhat exaggerated. Also shown in Fig. 3 is the total 
spectrum if the entire disc is assumed to emit an effective 
blackbody spectrum locally. The actual spectrum is seen to 
be much harder than that of a blackbody disc, with a greater 
spectral flux for e ^ 3kcV. 

Fig. 4 shows the summed disc spectra in hard X-rays 
that we calculate for models with accretion rates correspond- 
ing to / = 0.1, 0.2, 0.3, and 0.4. We flnd that each total 
photon spectrum can be flt extremely well by a power law 
spectrum with an exponential cutoff. 
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Figure 4. Summed spectra (solid curves) for accretion discs with 
/ = 0.1 (lowest), 0.2, 0.3, and 0.4 (highest). Also shown (dotted 
curves) are fits made over the 2-20 keV band assuming cutofi' 
power-law spectra. 



Table 1. Cutofi' power-law fits to summed disc spectra. 



/ 


eo (keV) 


r 


0.1 


1.23 


-0.19 


0.2 


1.48 


-0.09 


0.3 


1.75 


0.05 


0.4 


2.04 


0.18 



— oc (e/eo)" exp(-e/eo). (13) 

The values of the parameters eo and F for the flts, made 
over the 2-20 keV band, are listed in Table 1. The values 
of r are small, so the fits are close to simple exponential 
spectra. (Since we are fitting photon fiuxes rather than en- 
ergy fluxes, an exponential spectrum does not correspond 
to a single-temperature bremsstrahlung spectrum.) At low 
accretion rates, the F values are negative, so the power-law 
factors rise with photon energy. As shown in Fig. 4, the flts 
are accurate descriptions of the calculated disc spectra over 
the decade 1.5 ^ e ^ 15 keV in photon energy. 

3.3 Blurred disc spectra 

The simple summed spectra shown in Figs. 3 and 4 arc not 
completely accurate descriptions of the spectra observed by 
a distant observer. Even for a face-on disc, gravitational red- 
shift and the transverse Doppler effect, which depend on 
the radius r, will modify the observed spectrum, making it 
softer. For an inclined disc, the flux emitted at a given lo- 
cal photon energy at radius r will be spread over a range 
of observed energies by the Doppler effect. We take these 
effects into account to lowest order by using the method de- 
scribed by Chen, Halpcrn & Filippenko (1989) without the 
correction for light bending. 

The resulting "blurred" disc spectrum for an accretion 
rate corresponding to / = 0.3 and an angle of inclination 
i = 30° is shown in Fig. 5. Small wiggles in the spectrum 
are merely artifacts of interpolation between computed mod- 
els. We flnd again that the observed hard X-ray spectrum is 
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Figure 5. Observed hard X-ray spectrum (solid curve) for an ac- 
cretion disc witli / = 0.3 inclined at i = 30°, taking into account 
Doppler blurring and gravitational redshift. Dotted curve is a 2- 
20 keV fit assuming a cutoff power-law spectrum. Dashed curve 
shows the spectrum that would be observed for a blackbody disc. 



Table 2. Cutoff power-law fits to blurred disc spectra. 



/ 


i 


£0 (keV) 


r 


0.1 


10° 


0.95 


-0.45 




30° 


1.07 


-0.17 




50° 


1.22 


0.02 


0.2 


10° 


1.15 


-0.25 




30° 


1.29 


-0.05 




50° 


1.48 


0.15 


0.3 


10° 


1.37 


-0.06 




30° 


1.51 


0.09 




50° 


1.74 


0.26 


0.4 


10° 


1.61 


0.12 




30° 


1.76 


0.23 




50° 


2.03 


0.37 



all of these spectra, cutoff power laws provide excellent fits 
over approximately a decade of hard X-ray energies. The 
parameters for the 2-20 keV fits are listed in Table 2. Since 
— 1/2 ^ r ^ 1/2, all fitting functions arc dominated by the 
exponential term. For a given accretion rate, both eo and 
r rise monotonically with increasing inclination. Similarly, 
both parameters rise monotonically with increasing accre- 
tion rate for a fixed angle of inclination. 
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Figure 6. Observed hard X-ray spectra for accretion discs with 
/ = 0.1 (lowest), 0.2, 0.3, and 0.4 (highest) inclined at i = 10° 
(solid curves), 30° (dashed curves), and 50° (dotted curves). 

well fit by a cutoff powor-law spectrum. While F is nearly 
unchanged, the cutoff energy is reduced to eo = 1.51 keV 
since the spectrum is softer than the simple face-on summed 
spectrum. As shown in Fig. 5, the cutoff power-law fit is ac- 
curate over the decade of hard X-ray energies responsible 
for most of the luminosity of tlic amotion disc. For com- 
parison. Fig. 5 also sliows tlic blurred spectrum tliat would 
result if each point on the surface of tlie accretion disc emit- 
ted its effective blackbody spectrum. The observed spectrum 
is considerably harder than such a blackbody disc spectrum. 

Blurred hard X-ray spectra for accretion rates corre- 
sponding to / = 0.1, 0.2, 0.3, and 0.4 are shown in Fig. 6 for 
inclinations of 10°, 30°, and 50°. For a given accretion rate, 
the lowest inclination gives the strongest flux at low pho- 
ton energies, but at high energies Doppler blurring causes 
the strongest fiux to occur at the highest inclination. For 



4 DISCUSSION 

We see that tlie observed hard X-ray photon spectrum for 
a luminous LMXB accretion disc should resemble a soft 
(1 eo ~ 2keV) exponential decay multiplied by a shallow 
( — 1/2 ^ 1/2) power law. If the magnetosphere disrupts 
the disc before it reaches 3i?s, then the spectrum of the 
disc will be even softer. Of course, the accretion disc spec- 
trum only represents one contribution to the total observed 
spectrum of a LMXB. Other contributions may come from 
a boundary layer at the neutron star's surface or magneto- 
sphere and/or from a hot corona above the accretion disc 
(e.g., see Lamb 1991). 

Fitting EXOSAT spectra of five bright low-mass X-ray 
binaries with a cutoff power law plus a blackbody spec- 
trum, White et al. (1986) found higher values (0.6 ^ F ^ 1.6 
and 3 ;^ eo ^ 6keV) for the power-law parameters. Tliey as- 
sumed that the cutoff power law and blackbody represented 
the accretion disc and boundary layer, respectively. The dis- 
crepancy between their results and our models may in part 
be due to a single-temperature blackbody spectrum being a 
poor representation of the boundary layer and other emis- 
sion regions of LMXB. The difference between observed to- 
tal spectra and the accretion-disc spectra that we have cal- 
culated should give important spectral information about 
tliese poorly understood emitting regions. Further informa- 
tion may come from the time variability of the different com- 
ponents. 
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